Preequilibrium Emission of Light Fragments in Spallation Reactions 
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The ability to describe production of light fragments (LF) is important for many applications, 
such as cosmic-ray-induced single event upsets (SEUs), radiation protection, and cancer therapy 
with proton and heavy-ion beams. The cascade-exciton model (CEM) and the Los Alamos version 
of the quark-gluon string model (LAQGSM) event generators in the LANL transport code MCNP6, 
describe quite well the spectra of fragments with sizes up to ''He across a broad range of target 
masses and incident energies (up to ^ 5 GeV for CEM and up to ~ 1 TeV/A for LAQGSM). 
However, they do not predict the high-energy tails of LF spectra heavier than ''He well. Most LF 
with energies above several tens of MeV are emitted during the precompound stage of a reaction. 
The current versions of our event generators do not account for precompound emission of LF larger 
than *He. The aim of our work is to generalize the precompound model to include such processes, 
leading to increased predictive power of LF production. Extending the model in this way provides 
preliminary results that have much better agreement with experimental data. 



I. INTRODUCTION 

Emission of light fragments (LF) from nuclear reac- 
tions is an interesting open question. Different reaction 
mechanisms contribute to their production; the relative 
roles of each, and how they change with incident energy, 
mass number of the target, and the type and emission 
energy of the fragments is not completely understood. 

None of the available models are able to accurately 
predict emission of LF from arbitrary reactions. How- 
ever, the ability to describe production of LF (espe- 
cially at energies > 30 MeV) from many reactions is 
important for different applications, such as cosmic-ray- 
induced single event upsets (SEUs), radiation protection, 
and cancer therapy with proton and heavy-ion beams. 
The cascade-exciton model (CEM) Q version 03.03 
and the Los Alamos version of the quark-gluon string 
model (LAQGSM) 0, [| version 03.03 event generators 
in the Monte-Carlo n-particle transport code version 6 
(MCNP6) [3l , describe quite well the spectra of fragments 
with sizes up to *He across a broad range of target masses 
and incident energies (up to ~ 5 GeV for CEM and up 
to - 1 TeV/A for LAQGSM). However, they do not pre- 
dict well the high-energy tails of LF spectra heavier than 
"^He. Most LF with energies above several tens of MeV 
are emitted during the precompound stage of a reaction. 
The current versions of the CEM and LAQGSM event 
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generators do not account for precompound emission of 
these heavier LF. 

The aim of our work is to extend the precompound 
model in the codes to include such processes, leading 
to an increase of predictive power of LF-production in 
MCNP6. This entails upgrading the modified exciton 
model currently used at the preequilibrium stage in CEM 
and LAQGSM. It will also include expansion and exami- 
nation of the coalescence and Fermi break-up models used 
in the precompound stages of spallation reactions within 
CEM and LAQGSM. Extending our models to include 
emission of fragments heavier than ^He at the precom- 
pound stage already gives preliminary results with much 
better agreement with experimental data. 



II. THEORETICAL BACKGROUND 

Our models consider that a reaction begins with the In- 
traNuclear Cascade, referred to as the INC. The incident 
particle or nucleus (in the case of LAQGSM) enters the 
target nucleus and begins interacting with nucleons, scat- 
tering off them and also often creating new particles in 
the process. The incident particle and all newly created 
particles are followed until they cither escape from the 
nucleus, are absorbed, or, for nucleons, reach a thresh- 
old energy (roughly 10-30 MeV) and are then considered 
"absorbed" by the nucleus. 

The preequilibrium stage uses the modified exciton 
model (MEM) to determine emission of protons, neu- 
trons, and fragments up to *He from the residual nu- 
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cleus. We discuss the MEM in more detail below. This 
stage can have a highly excited residual nucleus undergo- 
ing dozens of exciton transitions and particle emissions. 
The preequilibrium stage ends when the residual nucleus 
is just as likely to have a An = -1-2 exciton transition as 
a Ati = — 2 exciton transition. 

In the evaporation stage, neutrons and protons in the 
residual nucleus can "evaporate," either singly or as frag- 
ments. The CEM evaporation stage is modeled with 
a modification of the Furihata's generalized evaporation 
model code (GEM2) Q, and includes light fragments up 
to 28Mg. 

During and after evaporation, the code looks to see if 
we have an isotope that has Z > 65 and is fissionable. 
If it is, and there is fission, then the code also allows 
evaporation from the fission fragments. 

There are two models that are not directly parts of 
this progression: coalescence and Fermi break-up. The 
INC stage only emits neutrons, protons, and pions (and 
other particles, when using LAQGSM at high energies), 
so the coalescence model "coalesces" some of the INC 
neutrons and protons into larger fragments, by compar- 
ing their momenta. If their momenta are similar enough 
then they coalesce. The current coalescence model can 
only coalesce up to ''He fragments, the same as the pree- 
quilibrium stage. The Fermi break-up is an oversimpli- 
fied multifragmentation model that is fast and accurate 
for small atomic numbers; in the current CEM model it 
is used when any residual nucleus or fragment has a mass 
number less than 13. 

The MEM used by CEM and LAQGSM [iHI calculates 
Tj , the emission width (or probability of emitting particle 
fragment j) as 

r,{p,h,E)= [ \i{p,h,E,T)dT, (1) 

where the partial transmission probabilities, A^, are equal 
to 

Xi{p,h,E,T) = ?f^^^.5R(p, h) X 

^ 

Eq. (2) describes the emission of neutrons and protons. 
For complex particles, the nuclear level density tu be- 
comes more complicated and an extra phase-space factor 
7j must be introduced: 

^^.«p3(|)p,-l. (3) 

Eq. (3) for 7^ is a only a rough estimation that we 
improve by parameterizing it over a mesh of residual nu- 
clei energy and mass numbers in our codes [13. As the 
MEM uses a Monte-Carlo technique to solve the master 
equations describing the behavior of the nucleus at the 
preequilibrium stage (see details in [l[), it is very easy 
to extend the number of types of possible LF that can 



be emitted during this stage. We generalize the MEM 
to consider the possiblity of emission up to 66 types of 
nucleons and LF, up to ^^Mg. As a starting point, for 
the inverse cross sections. Coulomb barriers, and binding 
energies of all LF we use the approximations adopted by 
GEM2 0. 

III. RESULTS 

Extending the Fermi break-up model to include heavier 
LF (up to j4 = 16) allows us to use it for nuclei with 
A > 12, yielding increased accuracy for reactions with 
light targets. Below are examples of calculations by CEM 
(with the expanded Fermi breap-up model) compared to 
experimental data d, 0| • 
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FIG. 1. Comparison of CEM03.03F results(solid red lines; 
"F" stands for expanding the range of LF in CEM) with ex- 
perimental data by Hagiwara et al. Q (open symbols) for 70 
MeV protons on a natural carbon target. We calculate for 
^^C only. 
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FIG. 2. Comparison of CEM03.03F results (solid red lines) 
with experimental data by Uozumi et al. @] (open symbols) 
for 190 MeV protons on a ^Be target. 

Results from the extended Fermi break-up model 
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achieve good agreement with experimental results for 
these light targets. 



200 MeV p + "AI -> ^^Li + ... 




T(Mev) 

FIG. 3. Comparison of CEM03.03F results (solid red lines) 
with experimental data by Machner et al. 01 (open symbols). 



1200 MeV p + "='Ni ^ 'Li + ... 
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1200 MeV p + "='Ni ^ 'Be + ... 
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FIG. 4. Comparison of CEM03.03F results (solid red lines) 
withd experimental data by Budzanowski et al. [1] (open 
symbols). 



Expanding the MEM to include heavier LF (up to 
^®Mg) yields increased accuracy for several reactions we 
investigate. Figure 3 compares new simulations from our 
expanded model with data by Machner et al. [sj for 200 
MeV p + 27 Al. 

We also find that the integral spectra for n, p, d, t, 
^He, and ^He (not shown), are not significantly impacted 
by this LF emission expansion. 

Figure 4 compares our results for 1200 MeV p + ""'Ni 
with new data by Budzanowski et al. 

Similar results for different LF spectra are obtained for 
several other reactions (see, e.g., |ll|). 

Our results indicate that expanding the MEM to in- 
clude LF preequilibrium emission significantly increases 
accuracy of the high-energy spectra compared to experi- 
mental data. 



IV. CONCLUSIONS 

Extending the CEM model to include emission of light 
fragments (LF) heavier than ''He (up to ^^Mg) in the 
preequilibrium stage results in significantly improved ac- 
curacy compared to experimental data for several reac- 
tions, especially in the high-energy tails of the spectra. 

Future work includes finding a global parametrization 
for 7^ , incorporating the expanded event generators into 
MCNP6, adding coalescence of heavier fragments, further 
exploring Fermi break-up, and upgrading the evaporation 
model. 

This study was carried out under the auspices of the 
National Nuclear Security Administration of the U.S. De- 
partment of Energy at Los Alamos National Laboratory 
under Contract No. DE-AC52-06NA253996. 
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